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Abstract 



We briefly summarize recent efforts in Garching for modeling stellar core collapse and post- 
bounce evolution in one and two dimensions. The transport of neutrinos of all flavors is treated 
by iteratively solving the coupled system of frequency-dependent moment equations together 
with a model Boltzmann equation which provides the closure. A variety of progenitor stars, 
different nuclear equations of state, stellar rotation, and global asymmetries due to large-mode 
hydrodynamic instabilities have been investigated to ascertain the road to finally successful, 
5o ' convectively supported neutrino-driven explosions. 
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Our neutrino-hydrodynamic simulations were performed with a conservative, high-resolution 
shock-capturing scheme that employs an exact Riemann solver and integrates the hydrody- 
namics equations with third order accuracy in space and second order accuracy in time. This 
code is coupled to a neutrino transport scheme that solves the energy-dependent moment 



equations of neutrino number, energy, and momentum of order (v/c) by making use of a 
variable Eddington factor closure which is calculated from a model Boltzmann equation [1]. 
In addition to radial derivatives, the moment equations are supplemented by (i^-dependent) 
terms that correspond to the advection of neutrinos with fluid motions in lateral direction 
(with vq being the lateral component of the velocity vector in spherical coordinates), and 
the equation of motion of the stellar plasma contains also lateral gradients of the neutrino 
pressure [2]. General relativistic effects are taken into account approximately by (i) a gener- 
alized gravitational potential (the radial part of which is constructed by a comparison with 
the Tolman-Oppenheimer-Volkoff equation, including terms due to fluid motion and neutrino 
effects) and (ii) relativistic redshift terms in the transport equations [1]. A comparison with a 
fully relativistic treatment in spherical symmetry yields very satisfactory results [3] . We have 
continuously modernized our treatment of neutrino interactions in the supernova medium, 
now including a variety of reactions and improvements (e.g., nucleon-nucleon bremsstrahlung; 
interactions between neutrinos of different flavors; correlations, weak magnetism and recoil in 
neutrino- nucleon interactions; electron captures on heavy nuclei according to recent calcula- 
tions by Langanke and coworkers) in extension to the standard treatment of neutrino-matter 
interactions in supernova simulations (see Refs. [1, 2, 4] and the references to a long list of 
original papers therein). 
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Figure 1: Left: Composition in the highly degenerate 1.38 Mq core of an 8-10 Mq progenitor [6] 
with O, Ne and Mg in the central region enclosed by an C-0 shell. (X) denotes the mass fraction of 
a representative neutron-rich heavy nucleus that appears in nuclear statistical equilibrium, and Y e is 
the electron fraction. Right: Mass trajectories, shock position (thick solid line that rises to the upper 
right corner of the plot, where it reaches the surface of the C-0 shell), and neutrinospheres (v e : solid; 
v e : dash-dotted; v^, u T) D T : dashed), and gain radius (dotted) for the collapsing O-Ne-Mg core. 
The mass trajectories are equidistantly spaced with intervals of 0.01 Mq. The (blue) hatched region 
is characterized by a dominant mass fraction of He. (The plots were taken from Ref. [7].) 

2 Results 

In the past two years we have applied our new neutrino-hydrodynamics code to core-collapse 
simulations in one and two dimensions for a broad variety of progenitors, spanning the range 
from O-Ne-Mg cores on the low-mass end (~8— 10 Mq) to 25 Mq on the high-mass side. We 
have investigated different input for neutrino-matter interactions and have recently also tested 
the effects of different nuclear equations of state. This work was mostly done in course of 
the PhD Thesis of R. Buras and the Diploma Theses of F.S. Kitaura Joyanes and A. Marek. 
Below some main results are briefly summarized. 

2.1 Explosion of an O-Ne-Mg core 

The main improvement of our new simulations of O-Ne-Mg core collapse — which we did 
so far only in spherical symmetry — compared to previous approaches is the more accurate 
treatment of neutrino transport and neutrino-matter interactions. Using the nuclear equation 
of state (EoS) of Ref. [5], we could not confirm [7] that a prompt explosion occurs as found 
in calculations with simpler neutrino treatment but also with a different EoS [8]. The shock 
is created at a mass coordinate of ~0.45 Mq (defined by the moment when the postshock 
entropy first exceeds 3 fee per nucleon) and stalls (defined by the time when the postshock 
velocity becomes negative) only 1.2ms later at ~0.8M Q , still well inside the neutrinosphere 
and before energy losses by the prompt u e burst could have contributed to its damping. We 
also do not find the powerful shock revival by neutrino heating as seen in Ref. [9] (cf. also 
Ref. [10]). Instead, the shock continuously expands due to the monotonically decreasing 
preshock density associated with the steeply dropping density at the interface between C-0 
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Figure 2: Top left: Density profiles of a sample of progenitor cores for stars with main sequence masses 
between 11.2 M© and 25 M according to Ref. [13] (13 M©, nl3Gio_ld), Ref. [14] (11.2, 15, 20 M©, 
sll.2Gio_ld, sl5Gio_ld, s20.0Gio_ld, respectively), Ref. [15] (a Type lb progenitor, slbGio_ld), and 
Ref. [16] (15 and 25 M©, 115Gio_ld and 125Gio_ld, respectively). The dashed lines show profiles of the 
latter models after we have evolved them to a central density similar to that of the other progenitors. 
Top right: Shock positions vs. time for the set of progenitors. No explosions were obtained in spherical 
symmetry. Bottom left: Mass accretion rates by the shock as functions of time. The steep density 
gradients at composition interfaces cause sudden drops of |M|. Bottom right: Luminosities of v e (solid 
lines), D e (dashed) and heavy-lepton z/s and P's (dash-dotted), measured by a comoving observer at 
500km, around the prompt v e burst (left panel) and during the post-bounce evolution (right panel). 

shell and He shell (Fig. 1). Towards the end of our simulation the mass accretion rate by 
the shock has correspondingly dropped to less than 0.03 M & s _1 . Although our simulation is 
not yet finally conclusive in this point, we see indications that a neutrino-driven wind begins 
to fill the volume between neutron star surface and shock and will lead to mass ejection 
with a rather low "explosion" energy and little Ni production, very similar to the outcome of 
simulations of accretion induced white dwarf collapse to neutron stars [11]. But the results 
of O-Ne-Mg core collapse can be sensitive to the properties of the nuclear EoS [12]. We plan 
to investigate this in future simulations. 
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Figure 3: Left: Adiabatic index T = (dlnP/dlnp) s for an entropy s = 1 fee per nucleon (and an 
electron fraction Y e = 0.4) for the equations of state of Ref. [8] ("Wolff", thin lines), Ref. [23] ("Shen", 
medium) and Ref. [5] ("L&S", thick). Middle: Shock positions and neutrinosphcric radii of v e as 
functions of time for collapse simulations of a 15 M Q progenitor (Model sl5a28 [20]) with the three 
nuclear equations of state. Right: Prompt v e burst (left panel) and post-bounce luminosities of v e 
(solid lines), v e (dotted) and heavydepton i/s and £/'s (dashed) for the core collapse simulations of the 
15 Mq star with the three different equations of state. (The plots were taken from Ref. [24].) 



2.2 One-dimensional core collapse for progenitors between 11 M Q and 25 M 

In spherical symmetry none of the investigated core collapse models, all evolved with the 
EoS of Lattimer and Swesty [5], developed an explosion until we stopped the simulations 
at typically 300 ms after bounce. This confirms the results of fully relativistic simulations 
by Liebendorfer et al. [17]. In spite of a rapidly dropping rate of mass accretion through 
the shock (Fig. 2) neutrino heating is not strong enough to reverse the postshock infall of 
the accreted matter, and the shock retreats in all cases, following the contraction of the 
neutrinosphere, after it had temporarily expanded to a maximum radius between 140 km and 
180 km (Fig. 2). The latter expansion was driven by the accumulation of accreted matter 
between shock and neutrinosphere during the early post-bounce phase when the accretion 
rate is very high. The bumps visible in the shock trajectories (Fig. 2) correspond to the 
density and entropy discontinuities at the interfaces between shells of different composition, 
where the shock experiences a sudden drop of the mass accretion rate. Note the astonishing 
similarity of the prompt v e burst, which is a consequence of the very similar structure of the 
inner ~1 M & of the Fe core for different progenitors before collapse (Fig. 2) and the fact that 
the structural similarity is still present at bounce. The large luminosity differences during 
the post-bounce evolution are caused by model-to-model variations of the structure near and 
beyond the outer edge of the iron core (cf. also [18]). First core-collapse simulations of 15 M Q 
and 25 Mq progenitors with electron captures on heavy nuclei being treated according to 
Ref. [19] suggest that these findings also apply when electron captures on nuclei instead of 
free protons dominate during the infall phase. 

2.3 Variations of the nuclear equation of state 

For a 15 M progenitor (Model sl5a28 [20]) we have investigated the influence of the nuclear 
EoS on core collapse, shock formation, and post-bounce evolution in three different cases. We 
ran simulations with 

(i) our standard EoS (Lattimer and Swesty, "L&S") [5], which is based on a compressible 
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liquid drop model and employs a Skyrme force for the nucleoli interaction (our choice of 
the compressibility modulus of bulk nuclear matter is 180 MeV, and the symmetry energy 
parameter 29.3 MeV, but the differences in the supernova evolution caused by other values of 
the compressibility were shown to be minor [21, 22]); 

(ii) a new relativistic mean field EoS (Shen et al., "Shen") [23] with a compressibility of 
nuclear matter of 281 MeV and a symmetry energy of 36.9 MeV; 

(iii) an EoS that was constructed by Hartree-Fock calculations with a Skyrme force for the 
nucleon-nucleon interaction (Wolff and Hillebrandt, "Wolff") [8] and has a compressibility of 
263 MeV and a symmetry energy of 32.9 MeV. 

The three equations of state show, for example, clear differences in the adiabatic index T = 
(d In P/d In p) s around and above the phase transition to homogeneous nuclear matter (Fig. 3). 

Significant differences were also found during the core-collapse calculations: 

(a) in the nuclear composition of the progenitor core and during infall; 

(b) in the maximum density reached at bounce, and the density in hydrostatic equilibrium 
afterwards; 

(c) in the shock formation point, which is shifted outward by ~0.05 M Q in case of the stiffer 
Wolff and Shen EoSs [23, 8]; 

(d) in the shock stagnation point (defined by the time when the postshock velocities start to 
become negative), which moves out by up to ~0.1M Q (or ~10km) for the stiffer EoSs; 

(e) concerning the existence and duration of a post-bounce phase where a shell with a sizable 
mass fraction of heavy nuclei occurs transiently outside of the nuclear core; 

(f ) in the radius of maximum shock expansion and subsequent shock contraction, which closely 
follows the contraction behavior of the nascent neutron star (Fig. 3) ; 

(g) in the peak luminosity during the prompt v e burst and the evolution of the post-bounce 
neutrino luminosities (Fig. 3); 

(h) despite of differences in details, an application of the Ledoux criterion suggests a qual- 
itatively similar behavior of convectively unstable regions in the neutron star and in the 
neutrino-heating layer. 

The publication of a detailed discussion of these results is in preparation (Marek et al. 2004). 
2.4 Two-dimensional models 

Convective overturn in the neutrino-heated postshock layer plays an important role during 
the supernova evolution and can provide the crucial push for delayed shock revival by an 
enhancement of the efficiency of neutrino energy transfer [25, 26, 27]. This is also seen in our 
recent simulations which employ a much improved treatment of neutrino transport compared 
to previous models. Although we found convective enhancement of shock expansion [28], we 
were, however, not able to confirm the successful explosions obtained in simulations with 
considerable approximations of the neutrino physics [29] . 

Rotation of the Fe core of the progenitor star, even at a moderate rate, supports strong 
postshock convection and brings the star much closer to an explosion by the neutrino-heating 
mechanism. In a simulation of a 15 M Q progenitor we assumed the Fe core to rotate rigidly 
with a rate of = 0.5 rad s _1 [28, 30], which is in the ballpark of results from stellar evolution 
calculations [31] and significantly slower than adopted in other recent works (e.g., by Burrows 
and collaborators [32] and Kotake et al. [33], see also their contributions to this meeting). 
Our rotating model showed shock expansion to a maximum radius of nearly 300 km along 
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Figure 4: Snapshots of the entropy distribution at t > 200 ms post bounce for two two-dimensional 
simulations with "Boltzmann neutrino transport" as decribed in Sect. 1. Yellow denotes highest 
entropies (around 15-25 /eb per nucleon), red, green, and blue 8 fee per nucleon) successively lower 
values. Left: Result at 245 ms p.b. for a 15Af Q progenitor whose Fe core was assumed to rotate 
rigidly with = 0.5rads _1 [28, 30], which is in the ballpark of predictions from stellar evolution 
models [31]. Equatorial symmetry was assumed. At the displayed time the shock reaches its largest 
expansion of nearly 300 km in the axial (in the plot vertical) direction, but lateron contracts again. 
Right: Situation in an 11.2 Mq star at 226 ms post bounce (when the simulation had to be stopped 
because of shortage of computer time) [36] . The shock has reached a maximum radius of 600 km with 
no sign of return. Global, violent low-mode bipolar oscillations of the postshock layer, which were 
possible in this simulation due to a 180° grid, support a weak explosion in this nonrotating model. 
Note the significant polar to equatorial deformation. 

the rotation axis (Fig. 4, left), whereas postshock convection in the same stellar model, but 
without rotation, was not strong enough to drive the shock much farther out than in spherical 
symmetry [28]. 

All 2D simulations in Ref. [28], however, were computed with a lateral wedge of about 
90°. In case of the rotating model, for example, equatorial symmetry was assumed. This 
restriction, however, constrains the possible modes of the flow pattern in the convective layer. 
The calculations in fact show growing power in large structures at later times after bounce, 
and models with a full 180° grid tend to develop a dominance of the I = 1, 2 modes in 
connection with global dipolar oscillations of the postshock layer, a phenomenon which is 
also seen in the three-dimensional case (cf. Refs. [34, 35] and the contribution by Scheck at 
this meeting). A first 2D simulation with a 180° grid and accurate neutrino transport for 
a (nonrotating) 11.2 Mq star reveals a dramatic difference compared to the 90° run of this 
progenitor in Ref. [28]: The shock expands continuously with no sign of return until the 
simulation was stopped at 226 ms post bounce. At that time the expanding shock is highly 
deformed and has reached a polar radius of more than 600km (Fig. 4, right). In the 90° 
simulation the shock was contracting again at that time [28]. 
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3 Conclusions 

The results reviewed above demonstrate the limits of current state-of-the-art supernova sim- 
ulations but also indicate the frontiers where ultimate success might be achieved in modeling 
convectively supported supernova explosions. Rotation of the stellar core, even at a moderate 
rate, has a strong supporting influence on shock expansion. The same is true for low-mode 
convection and global nonradial shock instabilities, which are captured by 2D simulations only 
in case of a full 180° grid and sufficient resolution to follow the amplification of perturbations 
between shock and neutron star by the "advective-acoustic cycle" [37]. The incompletely 
known nuclear equation of state adds a component of uncertainty to the problem. Our ID 
simulations with different EoSs show interesting and significant differences, but corresponding 
multi-dimensional simulations with accurate neutrino transport have yet to be done. 
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